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Near-IR line-strengths in elliptical galaxies: evidence for 
IMF variations? 
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ABSTRACT 

We present new relations between recently defined line-strength indices in the near-IR 
(CaT*, CaT, PaT, Mgl, and sTiO) and central velocity dispersion (ctq) for a sample 
of 35 early-type galaxies, showing evidence for significant anti-correlations between 
Can triplet indices (CaT* and CaT) and log do- These relations are interpreted in the 
light of our recent evolutionary synthesis model predictions, suggesting the existence 
of important Ca underabundances with respect to Fe and/or an increase of the dwarf 
to giant stars ratio along the mass sequence of elliptical galaxies. 
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1 INTRODUCTION 

During the last decade, the measurement and interpreta- 
tion of blue optical line-strength indices in the spectra of 
early-type galaxies in the field have revealed the existence 
of an apparent spread of mean ages (Gonzalez 1993; Faber 
et al. 1995; j0rgensen 1999) and element abundances ra- 
tios (Worthey 1998; Trager et al. 2000a), suggesting a va- 
riety of interpretations of scaling relations like the colour- 
magnitude or Mg2-cr relations (Bender, Burstein & Faber 
1993; Kuntschner 2000; Trager et al. 2000b; Vazdekis et al. 
2001). Since the picture from the blue is rather confused, if 
one wants to achieve a more complete understanding of the 
star formation history of these galaxies, it is necessary to 
look at other spectral regions in which the relative contri- 
bution of the distinct stellar types is very different. In this 
sense, the potential of the near-IR spectral range and, in 
particular, of the Can triplet is still almost unexploited. 

Since Ca is an a-element like Mg, it should be enhanced 
compared to Fe in giant ellipticals (Es). However, as sug- 
gested by several authors, Ca seems to follow Fe (O'Connell 
1976; Vazdekis et al. 1997; Worthey 1998; MoUa & Garcia- 
Vargas 2000; Vazdekis et al. 2001; Proctor & Sansom 2002). 
Even so, given that variations of the Fe line-strengths among 
Es are not negligible (e.g. Gorgas, Efstathiou & Aragon- 
Salamanca 1990; Gonzalez 1993; Davies, Sadler & Peletier 
1993; Kuntschner 2000), one should not expect the small 
variation of the Call triplet strength reported by previous 
work (Cohen 1979; Bica & AUoin 1987; Terlevich et al 1990; 
Houdashelt 1995). Furthermore, this result is difhcult to in- 



terpret in the light of previous stellar population models 
(Garcia- Vargas et al. 1998; Schiavon et al. 2000) which pre- 
dict a high sensitivity of the Ca 11 triplet to the metallicity of 
old, metal-rich stellar populations. Also, the absolute values 
of Can in Es differ from the model predictions (Peletier et 
al. 1999; MoUa & Garcia- Vargas 2000). 

With the aim of clarifying the above inconsistencies, 
during the last years we have developed a new stellar library 
in the near-IR spectral range (Cenarro et al. 2001a, hereafter 
CENOl) with a homogeneous set of revised atmospheric pa- 
rameters for the library stars (Cenarro et al. 2001b), deriv- 
ing empirical fitting functions that describe the behaviour 
of new line-strength indices for the Call triplet and the H 
Paschen series (CENOl; Cenarro et al. 2002) and other spec- 
tral features (Cenarro 2002, hereafter CEN02). Finally, in 
Vazdekis et al. (2002; hereafter VAZ02) we present a new 
evolutionary stellar population synthesis model which pre- 
dicts both the integrated indices and the spectral energy 
distribution for single stellar populations (SSPs) of several 
ages, metallicities and initial mass functions (IMFs). 

In this letter we present the first results for a spectro- 
scopic sample of 35 early-type galaxies. After a brief de- 
scription of observations and data reduction (Section 2), in 
Section 3 we describe the measurements of the new indices 
for the central regions of the galaxies and their relationship 
with the velocity dispersion. In Section 4 we discuss plausible 
interpretations of the data on the basis of new index-index 
diagrams derived from our model predictions. 
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2 OBSERVATIONS AND DATA REDUCTION 

Our sample consists of 35 carly-typc galaxies (E - SO) span- 
ning a wide range of absolute magnitudes (—22.5 < Mb < 
— 16.5 mag, using Hq = 75 kms~^ Mpc~^) and central veloc- 
ity dispersions (40 ao 370 km s"'^). Most of them are 
field Es, although a few galaxies from Virgo (9) and one cD 
in the Coma cluster arc also included. 

Long-slit spectroscopy was carried out during three 
nights in 1999 using ISIS at the 4.2 m WilUam Herschel 
Telescope (Observatorio del Roque de los Muchachos, La 
Palma), providing 2.9 A (FWHM) spectral resolution in the 
red arm (8355 - 9164 A). The slit (2arcsec width) was 
aligned with the major axis except for two SOs (along the mi- 
nor axis). Exposure times of 1200 - 2000 s per galaxy allowed 
us to obtain signal-to-noise ratios per angstrom from 43 to 
253 A^^ in the central spectra. We followed a typical spec- 
troscopic reduction procedure with REDmE (Cardiel 1999, 
see also CENOl), taking special care on the sky subtraction 
and the correction for fringing and telluric absorptions. The 
availability of error spectra for each galaxy frame allowed 
us to estimate reliable uncertainties in the measurements of 
the indices. The spectra were relative-flux calibrated using 
4 spectro-photometric standard stars (Oke 1990) observed 
several times at different air masses. Also, in order to cor- 
rect for small differences between the spectro-photometric 
systems of the galaxies and the model predictions, a sample 
of 49 stars (from B to late M spectral types) in common with 
CENOl were observed during twilights. They were also em- 
ployed as templates for velocity dispersion determinations. 



3 INDEX-logao RELATIONS 

In this section we present the behaviour of new line-strength 
indices in the ncar-IR spectral region (CaT*. CaT, PaT, 
sTiO and Mgl) as a function of the central velocity disper- 
sion ((7o) of the galaxy sample. 

CaT and PaT measure the strength of the Call triplet 
(AA8498, 8542, 8662 A) and three lines of the H Paschen se- 
ries. CaT* (= CaT - 0.93 PaT) is an index corrected for the 
contamination by the Paschen series in stars of the earliest 
spectral types (see CENOl for full details about their defini- 
tions). The sTiO index is a measurement of the slope of the 
continuum at the Ca II region. It is mainly governed by TiO 
absorption bands which are prominent in mid-late M types 
and observed in the integrated spectra of early-type gala^xies. 
Following CEN02, it is computed as the ratio between the 
CaT* pseudo-continuum values (C(A)) at the central wave- 
length of its reddest and bluest continuum bands, that is, 
sTiO = C(A8784.0)/C(A8479.0). Finally, Mgl measures the 
strength of the Mgl Une at A8807A (Cenarro et al. 2001c; 
CEN02). 

Central velocity dispersions for the galaxies were deter- 
mined using the MOVEL and OPTEMA algorithms (de- 
scribed in Gonzalez 1993) as explained in Pedraz et al. 
(2002). In each case, the template that was used was ob- 
tained as the mixture of 6 different spectral types (G5V, 
G7V, KG III, K2III, M3III and M5III) which minimizes 
intrinsic differences with the galaxy spectra. In order to 
avoid systematic differences between indices at diferent spec- 
tral resolutions, all the spectra were broadened up to the 
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Figure 1. Behaviour of the new indices with the central velocity 
dipersion {era) for the galay sample, at ct = 370kms^-i spectral 
resolution and corrected to the system defined by the models. Tlie 
lines represent error-weighted, least-squares fits of the form I = 
A + B log (TO to all data. Labels include the derived coefficients 
and their errors, the significance levels of Spearman rank-order 
tests (ar-s ) 1 the residual standard deviation of the fit (rms) and 
typical index errors for the galaxy sample (ctyp)- In the case of 
PaT, rms refers to the standard deviation w.r.t. the mean index 
value. 
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largest ao of the galaxy sample (~ 370 kms^^). The mea- 
sured indices correspond to a central aperture of radius 
Res/S (or larcsec for galaxies with Res < Sarcsec) and 
were corrected to the system defined by the models. See 
http : //www .ucm . es/inf o/Astrof /ellipt/CATRIPLET . html 
for a database with the indices and velocity dispersions. 

Figure 1 shows the obtained relations of the indices and 
the central velocity dispersions of the galaxies. Although, 
at first sight, the two Can triplet indices do not follow a 
neat linear behaviour with logtJo, the significance levels of 
Spearman rank-order tests (see the labels) show that nega- 
tive trends with velocity dipersion are significant. This re- 
sult is highly surprising since it is the first evidence for an 
anti-correlation between a metal-line index and tfie velocity 
dispersion. Note that classical metallicity indicators in the 
blue spectral range increase with ao (e.g. Mg2, <Fe>). In 
any case, it is worth noting that the spread of CaT* and 
CaT values is only ~ 5 per cent of the mean values. Prob- 
ably, this is the reason why previous works did not find 
significant variations of the Ca ll triplet in their galaxy sam- 
ples. For the sTiO and Mgl indices we find clear increasing 
trends with log (To, whereas we do not detect any significant 
trend for the PaT index. Only NGC4742 significantly de- 
parts from the mean PaT value, revealing a central young 
stellar population in agreement with other determinations 
in the optical (Gorgas et al. 1990; Trager 1997). 

Although a linear fit is representative of the general 
behaviour as a whole, different trends with velocity disper- 
sion are apparent. While the indices CaT*, CaT and sTiO 
of low mass Es (logtro ^ 2.20) are roughly independent of 
ao, galaxies with 2.30 ^ log ao ^ 2.50 depart from the above 
trend showing lower (CaT* and CaT) and larger (sTiO) val- 
ues. Also, some of the most massive Es (log ctq ^ 2.50) sig- 
nificantly deviate from the fit attaining the typical indices 
values of low mass Es. 



4 INTERPRETATION AND DISCUSSION 

To analyse the previous relations we make use of our 
SSP model predictions (VAZ02; CEN02), which were trans- 
formed to the spectral resolution of the data (370kms~^) 
using specific polynomials corresponding to their own age, 
metaUicity and IMF (see VAZ02). 

Given that the time evolution of the near-IR indices is 
virtually null for SSPs of all metallicities and ages ^ 3 Gyr 
(VAZ02, CEN02), we can consider that the age has a negligi- 
ble effect in this spectral range. Figure 2 shows the distribu- 
tion of the galaxy sample in the CaT*-sTiO and Mgl-sTiO 
planes, with symbols indicating different ranges of ao. The 
insensitivity of the indices to age is apparent in Figs. 2a,b. 
From Fig. 2a, one immediately can notice that, while low- 
mass Es (filled circles) can be roughly fitted by with SSPs 
of metallicity below solar (~ — 0.4dex), no age- metallicity 
combination can account for the low CaT* values of mas- 
sive Es (filled squares). The spread and location of the most 
massive Es (asterisks) will be further discussed. Moreover, 
typical random errors in CaT* arc not able to explain the 
spread among Es, suggesting that other parameters must 
be taken into account to explain the observed trend. In this 
sense, we consider the existence of i) non-solar abundance 



ratios (hereafter NSAR), and/or ii) systematic variations of 
the IMF, as possible solutions to the above inconsistency. 

NSAR could play an important role to interpret Fig. 2a. 
In fact, the high Mgi values in massive Es (Fig. 2b) con- 
firm the behaviour found for Mg in the optical. In this sce- 
nario, the position of massive Es in Fig. 2a implies that 
Ca should be underabundant in these galaxies, as was al- 
ready suggested by e.g. Vazdekis et al. (1997) using the 
Ca4227 Lick/IDS index. As discussed in VAZ02, the use of 
a-enhanced scaled-solar isochrones (Salasnich et al. 2000) 
predicts even larger CaT* values (by ~ 0.5A) for metal rich, 
old SSPs. Therefore, given that overabundances of other a- 
elomonts exist, and the contribution of Ca to the total metal- 
licity is lower than 0.5 percent, massive Es should have ex- 
tremely low abundances of Ca to account for the observed 
values. Isochrones with accurate degrees of enhancement 
for the different a-elements as well as stellar spectral li- 
braries with appropriate clement ratios are indeed necessary 
to tackle with the above problem. 

Could it be that massive Es are much more metal rich 
than the models, i.e., much more than two times solar metal- 
licity? In this case, the shift of the red giant branch to lower 
temperatures would predict a significant decrease in the 
CaT* and CaT values. However, one would expect then sTiO 
values much larger than those observed (and also higher Fe 
values in the optical), but the measured indices are in rear 
sonable agreement with the current model predictions. 

The alternative to interpret the data on the basis of the 
current model predictions is to assume that the IMF may be 
varying among Es. Since age can be considered a secondary 
parameter, in Fig. 2c we plot the SSPs model predictions 
for 17.78 Gyr (a different age docs not qualitatively alter 
the final conclusions) and different IMF slopes and metal- 
licities. In this new scenario, the data can be interpreted by 
means of IMF slope and metallicity variations, suggesting 
that low-mass Es (filled circles) exhibit lower metallicities 
and lower IMF slopes (lower dwarf/giants ratio) than mas- 
sive Es (filled squares). In the light of this interpretation, 
reading fj, and [Fe/H] values from the grid in Fig. 2c al- 
lows us to fit a one-parameter relation between ^, [Fe/H] 
and logCTo for Es in the sense that, the larger ao, the larger 
the metallicity and the IMF slope (Figure 3). Note, however, 
that very massive Es (asterisks) deviate from the above rela^ 
tion exhibiting lower metalhcities and slightly fiatter IMFs. 
To guide the eye, solid lines represent simple polynomial fits 
to the data on the /J [Fc/H] and [Fc/H] log cto planes, whilst 
the relation in the /i-log ao plane is readily derived from the 
two previous ones. Using Monte-Carlo simulations (follow- 
ing a procedure similar to that of Kuntschncr et al. 2001), 
we have checked that this relation is not driven by a com- 
bined effect of Poisson noise and non-orthogonal diagnostic 
diagrams, thus concluding that Es can indeed be explained 
using /Lt-[Fe/H]-log ao curves. 

Although the universality of the IMF is a highly con- 
troversial question (sec e.g. Eisenhauer 2001 for a review), 
there are theoretical arguments suggesting that the IMF of 
metal-rich star-forming regions must be biased towards low 
mass stars due to a more efficient cooling rate (e.g. Lar- 
son 1998). On this base, the universal IMF by Padoan et 
al. (1997) -whicli depends on the physical conditions of the 
star formation site- predicts such a behaviour. The obser- 
vational /^-[Fe/H] trend of Fig. 3c qualitatively agrees with 
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Figure 2. CaT*-sTiO and Mgl— sTiO diagrams for the galaxy sample corrected to the system defined by the models at 370 km s"-' 
spectral resolution. Panels (a) and (b) show SSPs model predictions with fixed Saljjctcr IMF slope (^t = 1.3). Age varies from 5.01 Gyr 
(dotted line) to 17.78 Gyr (dash-dotted line) with solid lines for intermediate values (A log [age(Gyr)] ~ 0.05). Metallicity spans from 
—0.68 to +0.20 as in the labels (dashed lines). Panel (c) shows SSPs model predictions at fixed age of 17.78 Gyr with varying power-like 
IMF slope {n = 0.3 — 3.3, see the labels) and metallicity (as in panel a). Different symbols indicate galaxies within distinct ranges of 
central velocity dispersion as it is shown in the key (b). Typical error bars for the whole sample are given. 
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Figure 3. The jU-[Fe/H]-log(To relation derived from Fig. 2c. 
Panels show the projection of the galaxies onto the planes ix— 
logCTQ (a), [Fe/H]-log (TO (b), and /i-[Fe/H] (c). Lines indicate 
least-squares polynomial fits to all data: (b) [Fc/H] = — 1.09 + 
0.46 logCTo; (c) ^l = 2.41 + 2.78 [Fe/H] - 3.79 [Fe/H]^; (a) is 
straightly derived from (b) and (c). Symbols are as in Fig. 2. 



these theoretical arguments, but it still needs to be related 
with the galactic mass (or logao). A time-extended (non- 
instantaneous) star formation history or several star-forming 
episodes could drive such a relation. If mergers and accretion 
of other galaxies with pre-enriched gas have been frequent 
phenomena during the evolution of massive Es, subsequent 
star formation episodes should increase the metallicity and, 
as a consequence, the IMF would gradually tend to produce 
a larger relative number of low mass stars. In this case, the 



infall of non-primordial, previously enriched, gas is needed 
to ensure the increasing metallicity. Although an extended 
star formation could be enough to explain the increasing 
metallicity, the addition of dissipativc hierarchical scenarios 
would emphasize such a behaviour. 

This interpretation implicitly characterizes the IMF as 
a time-dependent function (it gets steeper with time) in 
agreement with several previous work (e.g. Larson 1998), 
suggesting that the IMF in the early Universe could be bi- 
ased towards very massive stars. A prompt initial enrich- 
ment -straight consequence of a time-dependent IMF- could 
account for the observed a-cnhancements (Vazdckis ct al. 
1996). This could also explain several keys like the G-dwarf 
problem, the scarcity of Population III stars, and the high 
abundance of heavy elements in the intergalactic medium of 
clusters of galaxies (Larson 1998). In addition, the existence 
of a /i-[Fe/H]-log cto relation would also have important con- 
sequences for the interpretation of the blue spectra of Es. 
In particular, the <Fe>-log(To, H/3-logao, and the mass- 
luminosity (in B band) relations could only be reconciled 
by introducing an age sequence in the sense that the larger 
log (To the lower the mean age in the central parts (CEN02). 
Another problem could be visual-infrared colors. For exam- 
ple, for a model with p = 2.8, [Fe/H] = -1-0.2 and 12.6 Gyr 
one expects V - K = 3.52 (Blakeslee et al. 2001). This is 
just at the edge of the range of observational values for Es 
(e.g., V-K <. 3.50±0.05 for Es in Frogel ct al. 1978). How- 
ever, and again, a younger mean age for the central regions 
would help to match the observed values {V — K = 3.27 for 
a model with ^ = 2.8, [Fe/H] = +0.2 and 5.01 Gyr). 

The interpretation with a varying IMF actually revisits 
the classic debate about the existence of a dwarf-enriched 
population in the nuclei of Es, based on the strengths of the 
near-IR Nai doublet, the Call triplet and the FeH Wing- 
Ford band (Cohen 1978; Faber & French 1980; Carter et al. 
1986; AUoin & Bica 1989; Couture k. Hardy 1993). In partic- 
ular, the weakness of the FeH band was used as an argument 
against this possibility. Since these results are mainly based 
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on rather limited empirical synthesis models, they are not 
strong enough to directly exclude the possibility of a dwarf- 
heavy IMF. A proper calibration of its sensitivity to the 
stellar parameters (in particular to metallicity; see Carter 
et al. 1986) and its inclusion in modern stellar population 
models must be performed before extracting any conclusion. 
In any case, since several model uncertainties affect the ab- 
solute scale of the predicted index strengths (by ~ 0.5 A in 
the case of CaT*; see VAZ02), the derived ^i-[Fe/H]-log ctq 
relation must be considered on a relative basis. 

Finally, we speculate on the very massive Es labelled 
with asterisks in Figs. 2 and 3. Whereas all (except one) are 
boxy Es. with slow rotation and resolved cores, the rest of 
the sample arc mostly disky Es, fast rotators with power-law 
cores. According to Faber et al. (1997), disky Es are consis- 
tent with their formation in gas-rich mergers, whereas boxy 
Es could be the by-products of gas-free stellar mergers. In 
the context of our interpretation, this suggests that boxy 
Es should exhibit lower metallicitics (and, therefore, lower 
sTiO values) at the time that keep a more primordial, flat- 
ter slope IMF than disky Es (thus leading to larger values 
of CaT* and CaT), in agreement with what we observe in 
Figs. 1 and 3. Also, the fact that boxy Es seem to be older 
than disky Es (de Jong & Davies 1997; Ryden, Forbes & 
Terlevich 2001) favours the last interpretation. 

We want to note that, at the time of submission of this 
letter, Saglia et al. (2002) presented similar anti-correlations 
between Call indices and the velocity dispersion, using an- 
other large sample of high-quality spectra. Another paper by 
Falcon-Barroso et al. (2003) has been submitted with such 
anti-correlations for bulges of spiral galaxies. 

To conclude, more work in the areas of SN yields, stel- 
lar interior models and stellar libraries is needed to clarify 
whether it is possible to explain the current discrepancy be- 
tween Ca II measurements and models using non-solar abun- 
dance ratios. Until this is accomplished, and in the light 
of the present models, the only way out is to advocate for 
variations in the IMF. By no means we think that we are 
presenting strong evidence against its universality. Further- 
more, the scenario of a time-dependent IMF poses difficulties 
to explain other obsorvables. More modern, careful work us- 
ing indicators in the near-IR and the optical is needed to 
resolve the issue of a non-standard IMF. 
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